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FOREWORD 


Development of the ATLAS integrated structural analysis and 
design system was initiated by The Boeing Commercial Airplane Company in 
1969. Continued development efforts have resulted in the release and 
application of several extended versions of the system to aerospace and 
civilian structures. Those capabilities of the current ATLAS version 
developed under the NASA Langley Contract No. NAS1-12911 include the 
following: geometry control, thermal stress, fuel generation/management, 

payload management, loadability curve generation, flutter solution, 
residual flexibility, strength design of composites, thermal fully 
stressed design, and interactive graphics. The NASA monitor of this 
contract was G. L. Giles. The inertia loading capability was developed 
under the Army Contract No. DAAG46-75-C-0072. 

This document is one volume of a series of documents describing 
the ATLAS System. The remaining documents present details regarding the 
input data and program execution, data management, system design, the 
engineering method used by the computational modules, and 
system-demonstration problems. 

The key responsibilities for development of ATLAS have been 
within the Integrated Analysis/Design Systems Group of the Structures 
Research Unit of BCAC and the ATLAS System Group of the BCS Integrated 
Systems and Systems Technology Unit. R. E. Miller, Jr. was the Program 
Manager of ATLAS up to 1976 after which K. H. Dickenson assumed this 
position. The current ATLAS System is the result of the combined efforts 
of many Boeing engineering and programming personriel. Those who 
contributed directly to the current version of ATLAS are as follows: 
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ABSTRACT 


This document describes the automated design theory 
underlying the operation of the ATLAS Design Module. The 
methods, applications and limitations associated with the fully 
stressed design, the thermal fully stressed design and a regional 
optimization algorithm are presented herein. A discussion of the 
convergence characteristics of the fully stressed design is also 
included. Derivations and concepts specific to the ATLAS design 
theory are shown, while conventional terminology and established 
methods are identified by references. 
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NOTATION 


This list of notation complements the explanations provided 
in the text. Subscripted variables are only indicated with a dot 
in the following list and are explained in their proper context. 


A. 


b 

b 

B 

C 

Ck 

d 


Design variable; area 
Distance 

Design variation coefficient 
Design variation matrix 
Cosine of fiber direction (a) 
Value of equality constraint 
Value of inequality constraint 


D f D. Laminate constitutive matrix 

D 'i Lamina constitutive matrix 

Dn Design variation matrix 

e Subscript 

e • Exponent 

E,E. Moduli of elasticity 

E * Derivative of modulus of elasticity with respect to 

stress 

E. Element subsets (see ref. [10]) 

f , f . Applied stress 

f . Factor 

F, F. # F Allowable stress; (vector) 

F. Derivative of allowable stress with respect to thickness 

g^ Constraint function 


vii 



G r 6 . 


G,H 


i 

I 

j 

k 

*1 

l # m 


L,M 


MS, MS. 


n 

N 

N. 

o 


P. 


Q *ik 

R»R- 


R,R. 


R* 


S 




T 

Wi,W 

W 


Shear modulus 

Sum of inverted squares of allowable axial stress 

Subscript 

Unit matrix 

Subscript 

Subscript ; buckling coefficient 

Correction term 

Subscript 

Reciprocal of squared shear allowable stress 

Margin of safety 

Subscript 

Reciprocal of squared shear allowable stress 

End load 

Subscript 

Force 

Lamina compliance matrix 
Stress ratio 
Interaction function 

Derivative of interaction function with respect to 
thickness 

Sine of fiber direction (a) 

Thickness; design variable 
Allowable stress matrix 
Weight factor ; (vector) 

Weight 
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X 

y 

z 

a 

6,Y 

A 

V 

-X. 

A 

v 

P 


Subscript; vector 

Subscript 

Subscript 

Fiber direction; reciprocal of squared axial allowable 
stress 

Reciprocal of squared allowable stress 
Kronecker delta 
Difference operator 
Strain; (vector) 

Gradient 

Lagrange multipliers; (vector) 

Poisson" s ratio 
Density 

Stress; (vector) 
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1 . 0 INTRODUCTION 


The analysis segment of the structural design process has, 
for a reasonably long time, been supported by finite element based 
software. This support has resulted in practical and economical 
stress analysis of complex structures, like supersonic transport 
airplanes. As the design process, however, also requires 
repeated improvements and parametric changes of alternative 
structures, especially during the preliminary design phase, there 
is a great need for automated resize (design) methods. 

The ATLAS Design Module has been developed to satisfy these 
needs as far as strength and stability considerations are 
concerned. The methods have been developed and selected for 
application to large problems, (>1000 design variables). Some 
recent applications can be found in appendix C. 

It should be noted that the methods presented here are only 
intended as c omputation a l aid s in the design process. 

Engineering judgement is a required foundation for successful 
usage. This applies to the preparation of input data for the 
Design Module as well as interpretation and evaluation of the 
results. 

A judicious use, however, has been shown to produce both 
practical and economical benefits. It relieves the engineer of 
time consuming numerical tasks and makes him available for other 
design related activities (see refs. [ 1 ]-[ 3 ]) . 
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2.0 OVERVIEW OF THEOFY 


The general problem being addressed in this document involves 
the minimization of the structural weight (weight of the 
stiffness finite elements) , 


w = W(t x , . . . ,t n ) 


( 2 . 1 ) 


where t j represents H n M design variables subject to both the 
equality constraints 


f / . • • , t^) ~ C k (2.2) 


where k = 1 , . . . , rr<n, 

and the inequality constraints 


, . . . , fc^) ^ (2.3) 


The constraints are of stress type which involve strength and 
stability requirements that must be satisfied simultaneously. 

The problem identified above is clearly of optimization type. 
The solution of this problem will be approached from two 
viewpoints. One involves optimality criteria methods, whereas 
the other one concerns itself with ”math- programming” and local 
optimization. 
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2.1 OPTIMALITY VERSUS DIRECT SEARCH 

The optimality criteria and the direct search method can 
produce minimum weight designs. Caution is advised, however, as 
this is not always the case. 

The "math- programming '• methods (ref, [4]) are of direct 
search type which normally produce a "local" minimum. The 
computational state-of-the-art is, however, such that problems 
involving thousands of design variables are not practically 
solvable by these methods because of the cost and time factors. 
The "math- programming" methods are, therefore, limited to small 
and moderate size problems for which they work quite well. 

Recent developments involving the conjugate gradient method (ref. 
(5]) have shown promise for moderately large problems. 

The optimality criteria methods (ref, [6]) are of indirect 
type, i.e. no explicit minimization is required. Instead, a set 
of criteria is defined in such a manner that, when satisfied, a 
minimum or "almost" minimum is produced. Intuitive arguments are 
often used in selection of the criteria and consequently there 
will not be any guarantees that a "low-weight-design" has been 
achieved. Caution is clearly needed in interpretation of 
results. Experience shows, however, that optimality criteria 
quite successfully can be used to establish good resize 
procedures, and that for large problems these methods are the 
only feasible ones. 
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2.2 FULLY STRESSED DESIGN (F.S.D.) 

The "fully stressed design" method (F.S.D.) is of optimality 
criterion type. It represents the traditional way of sizing 
aerospace structures. That is to say, areas, thicknesses, etc. 
are chosen so that the applied loads yield stresses that are 
egual to the allowables. The redistribution of internal loads 
due to changes in structural properties is not considered until 
the next "stress-analysis," which consequently defines an 
iterative method based on repeated reanalyses . 

The method could be defined by the equations (2.2) when m=n. 
This yields a system of "n" equations in "n" unknowns. Solutions 
are, therefore, possible without any formal minimization steps. 
This is, however, in general a system of nonlinear equations for 
which convergence, uniqueness and a starting point must be 
evaluated separately for each application. 

The literature on automated design has an abundance of 
examples of the F.S.D. method producing undesirable results. The 
example illustrated herein is a square plate with a central 
circular hole (figs. 1-5) . This example is selected because the 
F.S.D. method is known to give poor results for this kind of 
problem. The first three figures show the change in thickness 
distribution. The fourth shows how the presence of several load 
cases tend to dampen the changes. Finally, figure 5 shows how 
stable the weight can be despite significant local changes. 

It is a well-known fact that the F.S.D. method works well for 
statically determinate structures. The main reason for this 
success is that a change in area or thickness for one member only 
changes the stress in that member and nowhere else. A matrix 
containing the derivatives of member stresses with respect to 
design variable would be diagonal for this case. 
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It is also true that for large structural models, like the 
ones used in the aerospace industry, the effect of changing a 
size variable is localized. This causes the matrix of the 
derivatives with respect to the size variables to be strongly 
diagonal . As shown in references [ 7 } and [ 8 ], this improves the 
chances of the F.S.D. to be optimum and to converge rapidly. 
Practical structures are quite often designed by multiple load 
conditions, a fact which also tends to improve the convergence 
behavior of the F.S.D. approach (see e.g. ref. [9]}. 

The F.S.D. method implemented in ATLAS is designed to be 
controlled by the user through input constraints and during 
execution by convergence specifications- The user can specify 
upper and lower bounds , input margins of safety, fixed values and 
regions to be excluded from resizing, (see ref. [TO}) and by so 
doing influence the convergence characteristics. 

It is finally noted that the "fully stressed design" method 
must be used with discretion and that results require judicious 
interpretation . When used in the preliminary design process it 
must be supported by experienced engineering personnel. However , 
it is still quite feasible to use it as a computational aid in 
the large problem environment. It is also quite frequent that 
the method produces consecutive designs with very minor weight 
changes even when member sizes exhibit significant changes (see 
figs. 4 and 5) . In conclusion, it is noted that the "fully 
stressed design" method is presently one of the few practical 
approaches to large problems. 
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Figure I . Thickness after One Cycle 

(Initial constant thickness) 
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Figure 2. Thickness after Five Cycles 

(Initial constant thickness) 
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Figure 3- Thickness after Ten Cycles 

(Initial constant thickness) 
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2.3 THERMAL FULLY STRESSED DESIGN (T.F.S.D.) 

The "fully stressed design" method is based on the assumption 
that internal element loads remain constant during resizing. A 
consequence of this assumption is that the following resize 
algorithm can be formulated 


A s= — a 
"new F "old 


( 2 . 4 ) 


where A new is the "new" value of the design variable, f is the 
applied stress, F is the allowable stress and A ^ the design 
variable value used when "f" was calculated. 

The algorithm shown by equation (2.4) has often been deemed 
adequate for resizing of structures subjected to mechanical 
loads. However, for structures with significant thermal 
influences it can lead to very slow convergence. The main reason 
for this is that stresses due to thermal effects quite often 
remain reasonably constant while the mechanical stresses change 
with changing size variables. An algorithm intended for 
situations involving significant thermal stresses has been 
developed and tested. Results and conclusions are shown in 
reference [11 J. A combination of thermal and mechanical stresses 
is considered and the resize algorithm is. 


new 


f 

- F-ff 


th 


old 


( 2 . 5 ) 


where the notation is the same as in equation (2.4) and "f t(l « 
represents the applied thermal stress. It is demonstrated in 
reference [11] how the use of this algorithm results in faster 
convergence than produced by the traditional "fully stressed 
design. " 
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2.4 REGIONAL OPTIMIZATION 


The successful application of “math-programming" methods in 
the large problem environment is contingent on the use of 
strategies that reduce a large problem to a set of smaller ones. 
That philosophy has been applied to the “composite optimization." 
This optimization is based on the method of feasible direction 
and is only used for structural models containing the composite 
elements C PLATE and C COVER (see sec. 4 and ref. [10]). 

Local optimization is used and total loads on the composite 
elements are assumed constant. The change in internal load 
distribution is, as for the F.S.D. method, introduced through the 
repetition of the stress analysis. 

The regional concept is introduced in such a manner that 
results, produced in the optimization, are valid for all elements 
in the region. These regions (subsets) are defined in the input 
data. 

The number of design variables available within the region is 
determined by the number of element variables. (All elements in 
a region have the same number of variables) . The number of 
constraints are, however, dependent on both the number of 
elements included in the subregion (see sec. 4) and the number of 
load cases involved. 

Criticality (relative value of constraints) of all the load 
cases is established during the initial phase of the 
optimization. This criticality is used to reduce the number of 
constraints, applied in the numerical search, to a number that is 
the same as the number of design variables. 

The problem defined in this manner is thus considerably 
reduced in size and guite feasibly solved by "math- programming" 
methods . 



3.0 OPTIMALITY CRITERIA 


The role of optimality criteria in automated design has been 
discussed extensively in the technical literature. It seems to 
be a common conclusion (see e.g. ref. [9]) that these kinds of 
methods, in the case of large structural systems, appear to be 
the most practical. Appendix C contains problems which presently 
cannot be solved directly by "math-programming" methods, and are 
consequently natural candidates for optimality criteria- based 
resizings. It is, however, also clear that for the size-range 
where "direct search" methods are applicable they constitute a 
quite natural complement to the optimality criteria approaches. 

The ATLAS Design Module contains two optimality criteria 
methods, the "fully stressed design" (F.S.D.) and the "thermal 
fully stressed design" (T.F.S.D.) • Both methods are based on 
repeated reanalysis, and convergence criteria can be of both 
weight and variable type. 

Both the F.S.D and the T.F.S.D. produce "margins of safety" 
that are used to update finite element properties. The resize 
algorithm used in both cases is. 


A 

new 


A old 

1+MS 


< 1+MS W 


(3.1) 


where "MS" is a calculated margin of safety and "MS| np " an input 
one. The design variables are denoted "A." It should be noted 
that all design variables are of "area" and "thickness" type. 

The details for each element are presented in appendix A. 

/The "margin of safety" defined here is, as the attentive 
reader already noticed, of an equivalent nature. It actually 
describes how a design variable should be changed so that the 
corresponding property value is just adequate for the applied 
loads. The more traditional definition of margin, on the other 
hand, results in a description of how the applied loads should be 
changed, so that applied and allowable stresses agree. 

The equivalent margins of safety, as calculated according to 
both the F.S.D. and T.F.S.D., are screened and the minima are 
used to establish an envelope of margins of safety. This 
envelope is used in equation (3,1). 
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3.1 FULLY STRESSED DESIGN THEORY 

The "fully stressed design" resize is based on three 
criteria. The Hill's criterion for strength (ref. £12]), an 
interaction criterion for panel buckling using input table 
allowables, and an interaction criterion for local buckling which 
uses calculated buckling allowables. 

The "fully stressed design" uses four kinds of allowable 
stresses. There are material allowables, average buckling 
allowables, gage-dependent table buckling allowables and local 
buckling allowables. The latter kind is calculated based on 
input "modulus tables" and "spacing" data. The buckling 
allowables are only applicable to plate- like elements. 

The calculation of "margins of safety" for buckling using 
"gage dependent allowables" involves two options. Both table 
allowables and local allowables are gage dependent and the 
margins can be calculated either as a first order approximation 
or by iteration to convergence. The change in allowable stress 
due to change in thickness is considered in both cases. 
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3.1,1 Str eng th Criterion 

Hill’s criterion (ref, [12]) is used to define failure for 
the ,, non-uniaxial H case. This criterion can be used for 
orthotropic materials and can, in its general form, be written 
as. 


F, VV 2LT yz + 2MT zx + 2 NT xy - 1 (3 ‘ 2) 


where '* 0." represents axial stress and " T.. '* shear stresses. The 
allowables-related variables are defined in the following manner. 



where ”F. M represents allowable axial stress. The shear related 
definitions are. 


2L 





(3.4) 


The variables "F.. M represent allowable shear stress. The 
plane case is traditionally established by an assumption of 
transverse isotropy; i.e. Fy=Fz • The resulting two-dimensional 
criterion becomes. 


33 



s 


(3.5) 




°x° y 

~2 

F 

X 


where F <1.0 when the criterion is satisfied. 


Extensive discussions of this type of criteria can be found 
in the literature (see e. g. refs. [12]-[16]). The advantage with 
this criterion is that it can be shown to reduce to the classical 
von Mises form for isotropic materials. It should be noted that 
only empirical arguments are available as far as explaining the 
analogy between yielding and failure. 
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3.1.2 Euc k linq Cri terion 

Buckling interaction criteria are used for the plate- like 
elements in the ATLAS system. Details can be found in appendix 
A. Three types of allowables are considered, average buckling 
allowables, gage- dependent table allowables and calculated local 
buckling allowables. A discussion of the former two is presented 
in reference [10] in the DESIGN input section. 

The local buckling allowables are of the form. 


E =k E (t)2 


(3.6) 


where "b" and the plate aspect ratios are defined in the "DETAIL 
DATA" in reference [10]. The assumption is that the geometry 
input defines rectangular plates on simple supports. The "k" 
values are calculated in accordance with the theory described in 
reference [17], (pages 348-406). 

The general interaction criterion is. 


R = 


(3.7) 


where fj to fj and e^ to ej are user specified parameters and 
"R. " represents ratios between applied and allowable stress. If 
the parameters ’are not input, the following default interaction 
expression is used. 



(3.8) 


This expression represents a parametric fit to results produced 
by the classical theory of buckling of simply supported, thin 
plates (see e.g. ref. [19]). A discussion of results using this 
type of interaction can be found in reference [18]. 
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3.1.3 Equivalent Mar g in of Safety 

It was shown in section 3.0 how the equivalent margin® of 
safety determine how the pertinent design variable should be 
changed in order to make the applied and allowable stresses 
equal. The general interaction expression can be written as. 


R o = R<t o> 


( 3 . 9 ) 


and the rrargin of safety is defined as. 


t + At 

o 


1+MS 


( 3 . 10 ) 


The intention is to change the design variable "t M so that. 


R = 1 = R(t Q + At) 


( 3 . 11 ) 


h Taylor expansion yields. 


1 = R (t Q + At) = R(t ) + At R' (t ) +o((At) 2 ) (3.12) 

h first order approximation and a substitution of the value for At 
yields , 


1 



MS 

1+MS 


t o R' 


( 3 . 13 ) 


and the equivalent margin of safety therefore- becomes. 
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MS 


( 3 . 14 ) 


1-R 


R-l-tR 8 

o o o 


where B* 0 is the first order derivative evaluated at t=t 0 . 


This is quite clearly an approximation that, for Converging 
or stabilizing design, becomes increasingly good. It is quite 
obvious that there are ranges for which equation (3.14) is 
meaningless. If that happens during the design, the module will 
revert back to the strength requirements and use those instead. 
Eguation (3.14) can be interpreted as a predictor that tries to 
correct for changes in, for example, allowable stress due to 
changes in plate thicknesses. 

Both of the interaction equations (3.7) and (3.8) can be 
written in the form. 


r = R(R X * W 


( 3 . 15 ) 


and derivatives of the type required in (3.14) have to include 
differentiation of R x , R y and R z with respect to "t" . The stress 
ratio R x is. 


R 


x 


N 


X ' 
tF x (t) 


( 3 . 16 ) 


where N x is endload and F x (t) represents the input table 
allowables. Differentiation yields. 



( 3 . 17 ) 


The variable F x ' symbolizes the derivative of the allowable 
stress with respect to thickness. These derivatives are 
obviously not available and the module, therefore, uses 
calculated "differences " to establish numerical "slopes” of the 
values in the tables. 
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When calculated local buckling allowables are used, the 
typical stress ratio K x becomes. 


R = 
x 


N b 
x 


t 3 E 


(3.18) 


The saire terminology as in equation (3,6) is used, and "N x " 
is the endload in the x-direction. The modulus of elasticity can 
be written as a function of applied «J) stress in the following 
manner. 


E = E(CT) 



(3.19) 


Differentiation with respect to the design variable (plate 
thickness) yields. 



(3.20) 


Finally, differentiation of 
”t M gives. 


equation (3.18) with respect to 



(3.21) 


This equation can be rewritten as 



( 3 . 22 ) 
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where E', as before, is the derivative of the modulus of 
elasticity with respect to applied stress. The value of this 
derivative is established from calculated differences of the data 
in the user supplied modulus tables. Allowable local shear 
stress and shear modulus tables are treated in an analogous 
manner. 

Expressions like equations (3. 17) and (3,22) must be 
evaluated in order to calculate the margin of safety (MS) in 
equation (3.14). This equation represents, as already discussed, 
a first order approximation. In order to improve on this 
equation, a local iteration is used as shown in section 3.1.4. 

The concept of "equivalent margin of safety" will finally be 
compared with the traditional approach. The interaction 
expression can again be written as. 


R — R(RxiRy/R X y) 


(3.23) 


which, with "increasing" load, becomes 


1 - R(R +AR ,R +AR ,R +AR ) 

x x y y xy xy 


(3.24) 


which up to a first order approximation is. 


1 = R + MS 


L 9R 

( x55 x 


» p 9& . p 9R 

y^R„ R xy8R 
J y xy 


) 


(3.25) 


If the "equivalent margin of safety" is used, the following 
expression is obtained. 


1 = R + 


/ MS e t VdH 
V 1+MS e ‘ I ^ 


(3.26) 
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where "MS e M represents “equivalent margin of safety. " The 
relation between the two kinds of margins becomes, from 
comparison of eguations (3.25) and (3. 26) , 


MS 


MS_ 

e 

1+MS 

e 


(1 + K r ) 


(3.27) 


where "Kj" is zero for the case of constant variables. For small 
values for MS e , 

MS ~ MS e (3.28) 
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3,1.4 Local Iteratio n 

De-termination of design variables, when "gage dependent 
allowables" are included, involves the so lution of the nonlinear 
equation, 

R(t) - 1 » 0. (3.29) 


This form immediately brings the "Newton-Raphson method" to 
mind (see ref. [ 20 ]) . The indicated method leads to a solution 
of the above equation by using the recurrence relation. 



R(t ± ) -1 
R' Wp 


(3.30) 


where "i" is the iterative index* A comparison shows that 
equation (3.30) is analogous to equation (3.14); or in other 
words the calculation of the "equivalent margin of safety" 
represents a first step in a "Newton-Paphson- like" iteration. 

Problems relating to divergence of the scheme suggested by 
equation (3.30) must be avoided. The convergence requirements 
will now be investigated by pursuing the analogy with "equivalent 
margins of safety." Equation (3.14) can be rewritten as. 


MS 


tR' 

R-l-tR' 


(3.31) 
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The general definition of the interaction equations yields 
the following relation, 

R'»RUt) < 0 (3.32) 

which results in the final form 

“S - - 1 + R-l^tR'l (3.33) 

The definition of margin of safety leads to the requirement, 

MS > -1 (3.34) 


which gives. 


R - 1 +|tR'| >0 (3.35) 


The default interaction is defined by equation (3. 8) which 
can be written as. 


R = R, + R 2 


(3.361 


Differentiation with respect to "t" yields 


I R't | > 4R 1 +R 2 > R (3.37) 


Substitution into (3.36) yields 

l't i K'(t|)J > * (3.38) 

This is a necessary requirement that is enforced during the 
iteration. The general interaction criterion in equation (3.7) 
leads to the same requirement shown by equation (3.38). 

Numerical experimentation has not produced a single case of 
divergence, once the said requirement is enforced. 
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3.2 THERMAL FOLLY STRESSED DESIGN THEORY 

The "thermal fully stressed design" algorithm is developed 
for structures subjected to both mechanical and thermal stresses 
at the same time. The thermal stresses are assumed to be 
significant compared to the mechanical ones, but are not allowed 
to exceed material allowable stresses. It is, in addition, 
assumed that mechanical internal loads and thermal stresses 
remain constant during resizing. It is, therefore, obvious that 
the thermal stresses cannot be allowed to be larger than the 
allowables, as that leads to an impossible situation. In a 
practical design situation one would expect that either a new 
material would have to be selected, or that a less severe thermal 
environment be found. 

The uni-axial case leads to the following expression. 


°n = 


a M A o 


n 


+ O- 
th 


(3.39) 


where 

®n ="new" total stress 

=mechanical stress (old) 
Oth “thermal stress 
A 0 = "old" area 
An^new" area 


which results in the "thermal fully stressed design" resize 
algorithm. 


M 


F '°th 


( 3 . 40 ) 


Here the total stress has been equated with the allowable. 
This expression illustrates quite clearly what happens when the 
thermal stress approaches the allowable. 


It is "tempting" at this point to extend the approach by 
using ratios between mechanical stresses and "modified 
allowables” in the interaction expressions used for the plane 
stress case. The analogy with equation (3.40) is intuitively 
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quite pleasing and has been explored extensively in the aerospace 
industry. However, a comparison between the traditional measure 
of criticality and the one indirectly suggested here, shows that 
in order for these to be equal the following equation must be 
satisfied. 


°M + °th 

F 


which is the case only when. 


a M “ F “°th (3 * 42) 


M 


F-a 


th 


( 3 . 41 ) 


As this obviously is true only when the ratios in equation 
(3.41) are equal to one, we find that the interaction expressions 
for “plate- like" elements require special attention. 

A comparison between the results produced by the T.F.S.D. 
method and the F.S.D. method are shown in figure 6. Different 
ratios for areas for different stress levels are plotted 
according to equation. 



V 


< F -°th ) ( °M +( W 


( 3 . 43 ) 


It is easily seen that quite drastic differences can be 
experienced. 
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Figure 6. Ratio of Areas Produced by the T.F.S.D. and FiS.D. Methods 
as a Function of Applied Mechanical Stress; CJ^ = „4F 


3.15 


3.2.1 Stre ng th Criterion 


The strength criterion used for the T.F.S.D. is the same as 
used for the F.S.D., (see eg. (3.5) ) , namely Hill*s criterion. 
If mechanical and thermal effects are kept separate, the 
criterion yields. 



The following simplified terminology has been used, 

N* = endload due to mechanical loads 
0 • = thermal stress 
F„ = allowable stress 

t = plate thickness (design variable) 

The resizing is now based on the requirement of zero margin 
of safety, (i.e. R=1.0). The following second order equation 

results. 



. 2 z , 

0(T„ + Yff xy - 1 



p 

— 


*1 

“j 

+ 2t 

“ N x°x 


< V 


+ 6 \ V™ x Ay 


+C® x 2 / 1 .M+0N y 2 + YN xy 2 . 0 
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The Greek letters are functions of allowables in the 
following manner. 



and the ’'new" thickness is calculated as the lowest positive root 
of equation (3.46). 


It is interesting to note how the traditional criticality 
approach compares to the one established for mechanical stresses 
and "modified allowables." The special case with one axial 
stress and shear will be investigated. The traditional way 
yields the ellipse. 


+ ( f xy +0 xy ) 2 _ 


p 2 

x 


2 
xy 


(3.46) 


and the "modified allowables" approach gives. 


£ ?y — = 1 (3.47) 

( f -<r 
\ xy xy ) 


A comparison between the two ellipses in figure 7 shows that 
the difference (the shaded region) changes with thermal stress 
and that quite unconservative designs can result. The "modified" 
curve consists of four different ellipses which, for small 
thermal stresses, would be a quite good approximation. The 
modified criticality ratio can be written as. 


f x 2 


( F x -°x) : 



(3.48) 
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Figure 7. Criticality Comparison; T.F.S.D. ( axial plus shear) 
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which, for small thermal stresses, becomes. 



( 3 . 49 ) 


It is furthermore expected that applied mechanical and 
allowable stresses will be reasonably close for this case which 
yields. 



( 3 . 50 ) . 


which could be interpreted as a good approximation. 
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4.0 COMPOSITE OPTIMIZATION 


The composite optimization in ATLAS is of a ‘‘local" type and 
is devised for the large problem environment* The method was 
developed for application to a composite SST with the number of 
design variables exceeding 6000. The composite membrane 
elements, CPLATE and CCOVER, are considered in the optimization. 
The CPLATE finite element has up to ten laminas of orthotropic 
material with user-defined fixed directions (see ref. [10]) • The 
CCCVEP is fcuilt up from two CPLATE elements. (Terminology can be 
found in appendix D) 

The structure to be optimized is considered to be divided 
into a number of regions (optimization problems) which are 
treated separately. These regions are defined by the input data 
and can be anything from one element to the whole structure. The 
design cf a structure consequently involves the repealed solution 
of weight optimization problems, each of which concerns itself 
wi+-h a small portion of the structure. 
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4.1 REGIONAL OPTIMIZATION METHOD 

The regional optimization method involves the solution of the 
following problem. For a given region (design subset of 
elements) El and an associated subregion (optimization subset) El. 
we have 



Figure 8. Regional Definition (E^ is a small part 
of the structure, and E^g can be several 
subregions) 

laminate strains and a set of initial lamina thicknesses. Based 
on this we want to find the lamina thicknesses M ti M # i=1 , . . . ,n 
where n<10, (analogously upper and lower surface for CCOVER) for 
all composite elements (of the same type and with the same number 
of laminas) in E ^ assuming regionally constant results. This is 
done in such a manner that minimum weight is produced without 
violation of any of the strength constraints in E^ . 
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4 . 1 « 1 C onst r aints 

The composite optimization is of strength type for which two 
optional types of constraints are considered; "Hill^s criterion'* 
and "the maximum strain criterion" (see refs. [ 14 ],[ 15] and 
[16]). 

The maximum strain criterion, which simply involves the 
comparison of applied strain to allowable strain for each 
component (2 axial and one shear) separately, uses the most 
severe one as the representation of criticality. The allowable 
strains are produced from input allowable stresses in the 
following manner. 




( 4 . 1 ) 


where "i" represents the lamina and "k" the laminate. The matrix 
» Q^" (in lamina reference frame) is defined as. 



1 

v 12 

E 11 

E 11 

v 12 

1 

E 11 

E 22 

0 

0 



(4.2) 


(see ref. [14], p 19). The allowable stress vector "F" is 
defined as. 


F xik 
F yik 
I F xyik 


( 4 . 3 ) 
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Tension allowable stress is used for both "F x jk " and M Fy|k " 
when tension allowable strain is to be produced, and analogously 
for compression allowable strain. 


The alternative to the maximum strain criterion is Hill’s 
criterion. This criterion has already been identified in 
eguation (3.5), but will be shown here in an alternative form. 
The criterion can be expressed as. 


<*ik T ik ff ik S X.O 


(4.4) 


where " <7jk i» represents the stress vector for the M i:th" lamina 
of the "kith" element. The allowables matrix "T^ " is defined 
as. 


T ik 


P xik 


2F 


xik 


2F 


xik 

0 


yik 


0 - 


(4.5) 


xyik 


The notation used here and in eguation (4.3) is 


F X jk = axial allowable stress in x-direction; 

F y [jf = axial allowable stress in y-direction; 

F xyik = shear allowable stress; 

in all cases for lamina "i" of element "k". Tension or 
compression allowable is selected in accordance with the sign of 
the comparable applied stresses. In order to use this criterion 
as a basis for resizing, we find that certain requirements with 
regard to allowable stresses must be satisfied. 
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The quadratic form in equation (4.4) must be positive 
definite. If that would not be the case, it would be possible to 
have stress fields for which the expression (4,4) remained 
negative, independent of how lamina thicknesses were changed. 

The determinant of the matrix , equation (4,5), is. 


det T ik = 


xyik 


xik 



( 4 . 6 ) 


It is easily seen from this expression that is not 
positive definite when Fyj^ >2F X || ( . In order to avoid this 
problem, it will be required that the larger axial allowable 
stress always defines the primary direction. 

The two criteria discussed here constitute the baseline for 
the constraints used in the optimization. The general form of 
the criteria for the maximum strain type is. 



e'ikj ; j=l, , . . ,3 | 

e*ikj Jj l 


( 4 . 7 ) 


ar.d for Hill’s, 


R ik£ ( a ik T ik°ik), 


( 4 . 8 ) 


where the index refers to load case and ”i” and "k" according 

to what was defined earlier. 

A detailed discussion of failure criteria for composites can 
be found ir references [12 ]-[ 16 ]. 
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4.1.2 Cr i t ic a lit y 

Each optimization problem is solved in an iterative manner. 
Each cycle in this iteration contains a screening phase and a 
solution phase. 

The screening phase consists of a search of subset E^. This 
search establishes the critical element and load case for each 
lamina "i”, identifying the following quantities. 


R. = max 
1 


where "ne" denotes the number of elements in E^ and M nl •’* is the 
number of load cases to be considered. The " s represent 

the value produced by the strength criterion used (Hill’s 
criterion” or "the maximum strain criterion”: see equations (4.7) 
and (4,8)). 

The screening is performed with the objective of establishing 
the strength constraints to be used during the solution phase. 
This procedure obviously requires all elements in Efcs to have the 
same number of variables (laminas) , and subsets containing CCOVER 
will be treated as two independent problems* 

The solution phase involves the minimization of weight with 
strength constraints defined by the screening. Thus, the defined 
optimization problem contains “n" design variables and ”n" 
constraints. Here ”n" is the number of laminas for any plate in 
E ks (upper or lower surface for CCOVER) . 

The screening is repeated after the local optimization is 
completed. If this does not result in "new" critical elements or 
load cases, the solution is complete, otherwise an additional 
cycle of optimization is performed and is repeated until the 
criticality order has stabilized. This is repeated a maximum of 
ten times for each problem. 


| R ikfc 


§ k— 1 f © • ,• f n6 7 J^**l f © « * § nl 




( 4 . 9 ) 
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4,2 LOCAL OPTIMIZATION 


The local optimization is of "math- programming” type that is 
concerned with the determination of a minimum weight design for a 
laminate. The lamina thicknesses (number of layers in each 
lairina) are the design variables, and there is one strength 
constraint present for each lamina. 

The starting point for the optimization is a set of laminate 
strains and an initial set of lamina thicknesses. The program 
establishes a starting point such that at least one lamina 
constraint is active. The starting point is determined from the 
information gained during the screening for criticality. The 
initial thicknesses sire determined from the input values, if 
Hill's criterion is used as, 

t is = V R imax' * t io ( 4 . 10 ) 


and if the maximum strain criterion is used as. 


t, - 
xs 


R imax 


xo 


( 4 . 11 ) 


where t.j- = starting point thickness, 
tj 0 = input thickness. 


and 


R -,*Tr,= „ - ma * 

imax ^ 


, x — 1 , « • . , n 


( 4 . 12 ) 


which establishes a starting point that is an acceptable design. 

The optimization is local in the sense that the total load on 
a laminate is assumed constant during the search. The 
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redistribution of internal loads due to changes in laminate 
stiffnesses is introduced through reanalysis. 
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4.2.1 Laminate Constitutive Relations 

A laminate, which in this context represents a OPIATE or an 
upper or lower surface of a CCQVER element, has in the general 
case anisotropic "membrane* properties. Each laminate is built- 
up of a maximum of ten directional laminas, each one with 
orthotropic elastic properties. 

For a laminate of total thickness »t" the following 
constitutive relation is applicable. 


a 


d e 


(4.13) 


where q represents gross "element" stresses and “c the gross 
"element" strains. (Element subscripts are dropped for 
convenience) . 

Each lamina has, in the local (lamina) reference frame, the 
following constitutive equation. 



* « 

D i c i 


(4.14) 


The theory is based on engineering strains and the strain 
transformation is, as easily can be verified in any elementary 
textbook. 



( 4 . 15 ) 
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where the transformation matrix "A " is. 



S 2 *■ SC 
c 2 -SC 


L-2SC 


2SC C 2 -S 2 J 


(4.16) 


Here "C" and "S” are# respectively, the cosine and sine of 
the angle (01, see appendix D) between the laminate reference 
direction and the lamina primary material direction. The 
transformation of stresses becomes consequently. 





(4.17) 


where " C\ n represents the stresses in lamina w i M referred to the 
laminate reference frame, and ” <7j” is referred to the lamina 
(local) frame. The total endload in the laminate, referred to 
the laminate frame can be written as, 

t i A?D^A i *e = tD ~t (4.18) 

i=-‘l 

which gives the following constitutive matrix for the laminate. 



n 


^2 tj a?d! a 4 


i x i a. 


i=l 


(4.19) 


where "C ” is the orthotropic constitutive matrix for lamina "i". 
The coefficients in this matrix are well-known and can, for 
example, be found in reference [14] page 18, or in appendix D. 
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4 . 2. 2 Lami na Stress and S t rain 

The input to the optimization procedure consists of a set of 
initial thicknesses and laminate strains. The optimization 
requires the repeated evaluation of stresses or strains as the 
design variables change. The lamina thicknesses "tj '* are 
perpetually compared to the initial values «t| 0 ", and stresses 
and strains are based on the assumption that total laminate load 
stays constant, 

N = tO = constant (4.20) 


This assumption provides the following expression for the 
laminate stresses. 





(4.21) 


where "r." indicates new values and w o" initial values. In 
addition, total strain compatibility is assumed in the laminate, 
yielding the lamina strains. 




A.D^D. 
in o 



(4.22) 


and the lamina stresses. 



d'.a.d' 1 d e\ 
i i n o o 


(4.23) 


Here only ’’tn’’ and ”E^ " are variables. It should be noted 
that, ir. the general case, a different element is ’’critical'* for 
each lamina (compare the screening phase in sec. 4.1.2). It is 
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quite clear from this situation that all the elements in the 
"design subset" Ej^must have the same number of laminas and 
identical orientations. 

In order to illustrate conditions surrounding a stable 
criticality situation, sufficiency in a special case will be 
shown. It is known that the initial lamina strains for two 
elements with the same orientations are related in the following 
manner. 


A i e ok > A i C ol 


( 4 . 24 ) 


A change in design variables results in the following 
difference between lamina strains. 


= € n VC<Wok t °fc- t ol D ol f ol ) 


If both elements have the same starting point. 


A i D ; lD ok (? ok- ? oi> 


( 4 . 25 ) 


which results in the same criticality as for the starting point. 
It is quite clear that this represents a very special case. 
However, it illustrates the phenomenon considered in the repeated 
optimization procedure, and indications are that initial 
criticality quite often is preserved. 
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4.2.3 Ma th-Progr am ming, M etho d 

The optimization is based on the method of feasible direction 
(see refs. [5] and [21]). The function to be minimized is the 
structural weight. 


n 

W • £ P i t i (4.26) 

i=l 


where , 

Pi = 

lamina 

density 


■t-i “ 

lamina 

thickness 


n .= 

number 

of laminas 


and the constraints are. 


9^ R^-l — 0 1 1—1 , . , * , n « 


(4.27) 


These constraints are of strength type, as discussed 
previously, and defined by equation (4.9) . 

The method of feasible direction (Zoutendi jk‘ s method 
described in ref. £4]) establishes a direction along which a step 
can be taken without violating the constraints, starting from a 
specific point in the design space. The feasible direction is, 
in this method, found by solving a linear programming problem in 
which the decrease in the structural weight "W" is maximized 
subject to constraints which insure feasibility; i.e. do not 
violate the constraints (4.27). 

The design variables "tj " are modified prior to the 
optimization, so that the largest constraint is equal to zero. 

The same normalization also takes place after the optimization is 
completed. Derivatives with respect to design variables, both 
for the objective function and the constraints, are established 
by using finite differences. 
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convergence of the optimization is considered to have been 
attained, if, in three consecutive iterations the relative and 
absolute change in the value of M W M is less than .001. The 
maximum number of iterations allowed is ten. 

Constraints can also be imposed directly on the design 
variables M tj ", Upper and lower bounds can be defined for each 
design variable, and the smallest lower bound is one layer 
thickness. No default exists for the upper bounds. 

It is, in addition, possible to define constraints that 
equate thicknesses for different laminas. These constraints are 
of the type. 


fc i t n y / * n l' n 2' ' n m 


where m<n-1. The resulting minimization problem involves less 
than "n" variables and M n H strength constraints of the type shown 
in equation (4.27). 
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APPENDIX A 


FINITE ELEMENT DESIGN 

This appendix contains a condensed description of each 
element type as seen from the standpoint of design. A detailed 
description of the finite elements is presented in reference 
[ 10 ]. 

Design variables, applied loads, allowables and criteria and 
algorithm information are presented for each element type. 
Additional information with regard to constraints (upper bounds, 
lower bounds, fixed data, etc.), allowables used, defaults and 
available options is described under "Design Input Data" in 
reference [10]. 
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ROD 


DESIGN VARIABLES: 

APPLIED LOAD: 

ALLOWABLES: 

DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 


One area at each end (2) 
Axial force 

Material allowable stresses 
None 

None 

F.S.D. and T.F.S.D. 


A„2 



BEAM 


DESIGN VARIABLES: 
APPLIED LOAD: 


ALLOWABLES: 

DESIGN PROPERTIES: 


REMAINING STIFFNESS 
PROPERTIES: 


DESIGN ALGORITHM: 


One AREA AND TWO SHEAR AREAS AT 
EACH END. (6) 

TWO BENDING MOMENTS AND ONE AXIAL 

Load at each end for the area. 

Two SHEAR FORCES AT EACH END/ ONE 
FOR EACH SHEAR AREA. 

Material allowable stresses 

Section moduli, two at each 
end; shear stress concentration 

FACTORS/ TWO AT EACH END. 

Area moments of inertia/ three at 
each end; changed in the same 
proportions as the area at the 

SAME END. 

F.S.D. 



SPAR 


DESIGN VARIABLES: 
APPLIED LOAD: 

ALLOWABLES: 

DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 


TWO AREAS FOR EACH CHORD PLUS 
WEB THICKNESS. (5) 

Average load for each chord, and 

EQUIVALENT SHEAR FLOW FOR SPAR-WEB. [£> 

Material allowable stresses/ plus 

SHEAR BUCKLING ALLOWABLES FOR THE 
WEB. 

Stiffener area and spacing 
Lumped cap areas/ changed by the 

SAME FACTOR AS THE WEB. 

F.S.D. and T.S.F.D. 



STIFFENER AREA CHECKED FOR SHEAR CARRIED AS DIAGONAL TENSION, 



PLATE 


DESIGN VARIABLES: 
APPLIED LOAD: 
ALLOWABLES: 


DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 

DESIGN INTERACTION: 


SIZING: 


ASSUMPTIONS: 


Plate thickness and two smeared 
STIFFENER THICKNESSES/ (3) 

Membrane endload/ two axial and 

ONE SHEAR COMPONENT, 

Material allowable stress 
Axial buckling allowable stress ) 
Shear buckling allowable stress j 
Local axial buckling stress ) 
Local shear buckling stress j 

2 STIFFENER RATIOS 

None 

F.S.D, and T.F.S.D, 

Hill's criterion for strength 
User defined interaction for 

BUCKLING, 

Predictor approximation/ 

Local iteration; both for gage 

DEPENDENT ALLOWABLES 

Tension is neglected in buckling 

INTERACTION - STIFFENERS USE THE 
SAME ALLOWABLES AS THE PLATE 


|£>Input gage dependent tables, Calculated, 



GPLATE 


DESIGN VARIABLES; 
APPLIED LOAD; 
ALLOWABLES; 


DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 

DESIGN INTERACTION: 

SIZING: 

ASSUMPTIONS: 


Up to five thicknesses? one at each 

CORNER AND ONE IN THE MIDDLE, (5) 

Membrane end loads only, two axial 

AND ONE SHEAR COMPONENT. 

Material allowable stress 
Axial buckling allowable stress) 
Shear buckling allowable stress) 
Local axial buckling stress ) 
Local shear buckling stress j 

None 



Bending thicknesses, unchanged. 

F.S.D. 

Hill's criterion for strength. 

User defined for buckling. 

Local iteration for gage dependent 
allowables 

Average membrane stresses are 
representative; bending stresses 

ARE SMALL AND NOT DESIGNING, 


> Input gage dependent tables. Calculated. Based on 

AVERAGE GAGE 



SPLATE 


DESIGN VARIABLE: 
APPLIED LOAD: 
ALLOWABLES: 

DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 

SIZING: 


One thickness (1) 

Equivalent shear flow 

Material allowable stress 
Shear buckling allowable stress 

None 

None 

F.S.D. and T.F.S.D. 

Predictor approx imi at ion 


|£>Input gage dependent tables. 



CP LATE 


DESIGN VARIABLES: 

APPLIED LOAD: 

ALLOWABLES! 

DESIGN PROPERTIES: 

REMAINING STIFFNESS 
PROPERTIES: 

DESIGN ALGORITHM: 

DESIGN CRITERIA: 

SIZING: 

ASSUMPTIONS: 


Number of layers for each lamina 
(max. 10) 

Laminate strain, three components 
Composite material allowables 
None 

None 

Composite optimization 
Hill's criterion or the maximum 

STRAIN CRITERION 

Math-programming; feasible direction 
Total load in laminate remains 

CONSTANT DURING RESIZING 
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REMAINING ELEMENTS 


COVER: 

Analogous to two PLATE: s. 

BRICK: 

Margin of safety calculated according 
to Hill's criterion 

SCALAR: 

Not designed 

SROD: 

The same as ROD 

CCOVER: 

Analogous to two CPLATEis 
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APPENDIX B 

CONVERGENCE OF THE FULLY STRESSED DESIGN 

Application of the "fully stressed design" method to 
practical structures normally leads to concern with regard to the 
convergence characteristics and "closeness" to minimum weight. 
Section 2.2 contains illustrations of the convergence behavior 
under adverse conditions. 

A discussion of the fully stressed design method is presented 
in references [7] and [8], Both papers discuss the problems 
relating to convergence. Some attention has been given to the 
question of verifying optimality for practical problems. No 
criterion has, however, been developed for applications to 
structures with many design variables, but some of the analytical 
aspects of the formulation are quite important even from a 
practical standpoint. 

The fully stressed design method is of iterative type and its 
application in a practical design situation requires 
determination of the following: 

•Is the method converging? 

• When should the iteration be terminated? 

• If it does not converge, what should be done from a 
constraint standpoint, so that the modified problem will 
converge? 

• How close to minimum weight is the solution? 



The fully stressed design, at the final point of convergence, 
satisfies the following set of equations. 




(B.l) 


where i=1, . . . ,n and “n" is the number of design variables 
present. The variables are. 


A j = design variable, 

Pj = equivalent design load, 

Fj = allowable stress. 


The general design problem, on the otherhand, is, 

n 


minimize 

subject to the constraints: 


»- E 

1=1 


w i A i 


(B.2) 


(B. 3) 


g i * A i"F7-° 

where i=1, n. Before the conditions surrounding optimality 

are described a few concepts are defined. 


The problem of convergence and optimality of the fully 
stressed design will be studied in a range where the changes in 
the design variables are reasonably small. A linear expansion of 
the critical forces yields. 


(k) _ p (k-1) 
i “ P i 



aP i (k-l) 




(B. 4) 


where the superscript denotes the iteration number. The design 
variation matrix, B=[ b jj ] is defined by differentiation of 
equation (E.1) yielding. 



9P i 

3aT 


(B . 5) 
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The design vector has the design variables as components 

and the superscript indicates the iteration number. The "design 
modification" vector is. 


A (k) = A (k+l)_ A (k) 


(B. 6) 


which also can be written as. 


AA* n> = B-B 0 . . .B AA (o) 
l z n 


(B. 7) 


where B n should be evaluated at a point between k-1 and k. The 
iterative process will converge if. 


D - B.B_ . . ,8_ -*• 0, when n -► oo (B.8) 

n 12 n 

A sufficient condition for convergence is that for the 
spectral norm of D n the following holds. 


|D n | -+ 0, when n-too (B.9) 

which is the case if, 

jB^j < 1, for 1=1, (B. 10) 


The spectral norm is defined as 



where |x| is the Euclidian norm of x. It is quite clear from 
(B. 11) that all the coefficients of Bj must be less than zero. 

The Euclidian norm for every column vector in Bj must also be 
less than zero. The relations described here do not constitute a 
practical method for determination of convergence requirements, 
but provide a background for tetter understanding of the 
convergence problems enountered with the fully stressed design 
method. 


The relation between the "fully stressed design" and the 
minimum weight design is of great concern. When these two 
designs coincide, they will both be solutions to equations (B.2) 
and (B. 3) . It is a well-known fact that for statically 
determinate structures equation (B.5) becomes. 



(B.12) 


in which case the right hand side of equation (B.3) becomes a 
constant and the minimization becomes a linear programming 
problem. The solution is always at the vertex of "n" of the 
"constraints". As there are only "n" constraints, we find that 
the fully stressed design represents a minimum. 

Eor indeterminate structures the "internal loads" are 
nonlinear functions of the design variables, and the minimization 
becomes a nonlinear programming problem. The solution is not 
necessarily a vertex, and the fully stressed design will 
consequently not always be a minimum. 

It is of great importance, in the practical situation to 
determine how close a fully stressed design is to the optimum. 

The Kuhn-Tucker optimality condition (ref. £5]) is. 


n . 

VW = 2 A j V9j ' (B. 13) 

3=1 

where y is the "gradient" symbol, and n is the number of active 
constraints. Equation (B. 13) must be satisfied at the optimum 
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point, which if coincident with the fully stressed design, would 
have to he a vertex. In that case. 




k=l 


hq. ’ 
dtl * 


W i 


(B. 14) 


and 


5 ^ ~ 0 for 1 / • ■ • / n 


(B. 15) 


Differentiation of the constraint equation yields. 





(B. 16) 


where 6jj is the Kronecker delta, and equation (B.14) can be 
written as. 


(i-b t ) A = W 


(B. 17) 


where I is a unit matrix, B T is the transpose of the matrix 
defined in equation (B.5). The variable X is a vector with the 
Lagrange multipliers as components and vf“ is a vector containing 
the right hand sides of all the equations similar to (B.14), The 
optimality condition is. 


A = (I-B T ) _1 w >0 (B. 18) 


which should be interpreted as every Component having to be 
positive. The requirement of positive Lagrangian multipliers is 
part of the Kuhn-Tucker conditions. 
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It is known that for a converging fully stressed design, the 
matrix I-B^ is nonsingular and the optimality condition can be 
written as, 

X* [i+B T +B T 2 + . ...]w> 0 (B. 19) 


The series inside the brackets converges rapidly when the 
coefficient b|j are much smaller than unity. Finally, if the 
matrix in brackets becomes diagonally dominant with positive 
diagonal values, the fully stressed design is optimum. 

Equation (B.18) shows that verification of the optimality of 
the fully stressed design requires computation of matrix B at the 
point of convergence. An approximation of this matrix can be 
produced through the use of finite differences. The coefficients 
would be calculated in the following manner. 


A (k+l)_ A <k) 

b • • ~ —3 3 

AA. 


(B. 20) 


This approach is presently the only one available, and it can 
be seen that it would be feasible only for very small problems, 
as up to "n" reanalyses could be required. 

The state-of-the-art is presently such that the convergence 
and optimality questions can he addressed only from inituitive 
standpoints. Experience seems to indicate that for rapidly 
converging fully stressed designs, the outcome roost likely is an 
optimum. For slowly converging designs, however, optimality must 
be investigated. 
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APPENDIX G ' 

SIGNIFICANT AUTOMATED DESIGN APPLICATIONS 

This appendix describes a few applications from the field of 
preliminary design. All examples have been analysed and designed 
using ATLAS. These structures are dominantly of SST-Iike 
character, with one exception. All belong to a size- range that 
presently renders math- programming global optimization methods 
impractical. 

representative size and computing requirements are shown in 
figure C. 3. A more detailed description of the three airplanes 
presented here can be found in references [ 1 ]-[ 3]. 

The examples are shown because of their sizes and because 
there are no implications with regard to convergence. The weight 
for each cycle is illustrated in figures C.1 and C.2 as a 
representation of what can be expected. 

Ihe two SST models, the Arrow Wing and the NST, as well as 
the "freighter” require further study in order to establish 
convergence and minimum weight. 
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Figure C.|. Arrow Wing Structural Weight 
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Figure C.2. NST Structural Weight ( with comparison 
to Multi- Mission Freighter) 
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IDENTIFICATION 

SIZE \ 

NO. 

NAME 

NUMBER 

DOF 

NUMBER 

D.V.. 

i 

METAL; 

ARROW WING 

9000 

20000 

2 

NST; base-line 

8000 

16000 

3 

NST; substructure 

8000 

10000 | 

4 

MULTI-MISSION 

FREIGHTER 

9000 

20000 

5 

COMPOSITE ; 

ARROW WING 

9000 

6000 


ORDER OF MAGNITUDE : Degrees of Freedom (DOF) 

Design Variables (D„V.) 


9000 

8000 

7000 

6000 

5000 

4000 

3000 

2000 

1000 


Design 

Stress analysis 


NO.. 1 2 3 4 5 

Figure C.3. Computer Time Requirements for One Cycle 
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APPENDIX D 

COMPOSITE TERMINOLOGY 

The ATLAS system has two composite elements that are 
identified by the names CPLATE and CCOVER. 

The CPLATE is a triangular or quadrilateral laminated 
constant-strain membrane plate. Each plate element (laminate) is 
comprised of up to 10 orthotropic laminas. Each lamina is 
defined by, 

• fiber direction (a) 

• material properties 

• temperature 

• number of layers 

The layer thickness is part of the composite material 
definition (see ref - [ 10 ]) . 

The CCOVER element is built up from two CPLATE elements 
separated by rigid posts. 
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*. Lnirn layers 


AREF-U. 

( or AREF-L ) 



1 (Fiber Direction) 

Reference 

Direction 


X U (° r *l) 


Figure D. I . Typical Lamina ** i M 


The orthotropic constitutive matrix for lamina T is* 


D S = 


E ll /( 1 -V12V21) 


Vl2E 2 2 / (l-Vi 2 V2l) 0 


V 1 2 E 22/< 1 ' V 1 2 V 21) E 2 2/a-V 1 2V2 1 ) 0 



0 


Q 



where engineering strains are used and where , 

V 21 e U = v 12 E 22 
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APPENDIX E 


DESIGN THEORY SUMMARY 

The ATLAS Design Module contains four different functions; 

• Fully stressed design; F.S. D. ; 

• Thermal fully stressed design, T.F.S.D.; 

• Composite optimization; 

• User specified changes. 

Each o£ v these automated design methods is of reanalysis type; 
i.e. changes in internal loads due to changes in thicknesses and 
areas are considered through a "new" stress analysis. The Design 
Module execution is consequently a part of the design process 
which, in addition, involves the stress analysis modules. It is 
primarily intended to be a part of the design process which 
involves the stress analysis modules. It is primarily intended 
for preliminary design applications. The first three functions 
above are all of iterative type. 

The "closeness" to minimum weight and the convergence 
characteristics are, of course, of great importance. It is, 
however, presently not possible to give general rules for the 
convergence behavior. The successful application of these 
methods will, therefore, ultimately have to depend on engineering 
judgement. 

The Design Module input data has been structured in such a 
manner that the convergence behavior of the structure can be 
changed during the design process. The input data contain the 
following constraints; 

• upper bounds 

• lower bounds 

• fixed data 

• input margins of safety 

• regional exclusion from resizing 

All these data sets can be modified or changed through an 
update capabi lity that allows for partial input data sets . This 
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feature, in addition to the fourth function specified above, 
allows the user to constrain the process or to change the 
direction in which a •'design” is developing. 

The F.S.D. and the T.F.S.D. methods produce equivalent 
margins of safety that are based on the resize algorithm. 


A new (1+MS) (1+MS inp ) 


(E.l) 


where MS is calculated and MSj np is an input margin of safety. 

The design variables involved are thicknesses and areas. 
Additioral finite element properties are changed in the same 
proportions as the associated design variables are modified. 

The F.S.D. method considers both strength and buckling 
requirements. Hill's criterion is used for strength, and user 
specified interaction criteria are applied for buckling. These 
two criteria are described in sections 3.1.1 and 3.1.2. Margins 
of safety produced by these criteria are based on the assumption 
of constant internal loads. The buckling criterion, however, can 
contain gage dependent allowable stresses. These can be of input 
type, table allowables, or they can be calculated local 
allowables. Eoth cases lead to a resize requirement that 
involves solutions of nonlinear equations. The principle of the 
solution of these equations shown in section 3.1.3. 

The T.F.S.D. method considers only strength requirements 
which are based on Hill's criterion. The method is intended for 
preliminary design of structures that are subjected to both 
mechanical and thermal loads at the same time. The method is 
described in section 3.2. It is based on the assumption that 
mechanical loads and thermal stresses remain unchanged during 
resizing. The thermal stresses are considered to be of the same 
order of magnitude as the mechanical stresses. 

The composite optimization is intended for finite element 
models containing CPLATE and CCOVER elements. This optimization 
is of local type, and lamina thicknesses are design variables. 
Regions containing the same type of elements with the same number 
of laminas and user-specified fixed fiber directions can be 
defined in the input data. Section 4.1 contains a description of 
these definitions and reference [10] describes the details. 
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The optimization is performed for minimum weight. The method 
of feasible direction is used and the constraints are of strength 
type. These strength constraints can be based on either Hill’s 
criterion or on the "maximum strain criterion." The local 
optimization is based on the assumption that total laminate loads 
remain constant. The changes in internal loads are introduced in 
the "next" stress analysis. 

The function that changes the finite element properties 
according to user specifications is intended as a means of 
introducing design requirements other than strength and stability 
based element criteria. This function consists of two 
algorithms, one that allows factoring of element subsets (flutter 
changes, e.g.) and one that allows individual changes to the 
finite element properties. This provides a capability for model 
changes that does not require "new" finite element input data. 
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